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Higher moments of event-by-event net-proton multiplicity distributions have been 
applied to search for the QCD critical point. Model results are used to provide a 



ly-j ■ baseline for this search. The measured moment products, ko 2 and So of net-proton 

distributions, which are directly connected to the thermodynamical baryon num- 
ber susceptibility ratio in Lattice QCD and Hadron Resonance Gas (HRG) model, 
■ are compared to the transport and thermal model results. We argue that a non- 

monotonic dependence of no 1 and So as a function of beam energy can be used to 
search for the QCD critical point. 
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^ ' I. INTRODUCTION 
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Lattice QCD calculations predict that a cross-over from the hadronic phase to the Quark 
Gluon Plasma (QGP) phase occurs above a critical temperature with zero baryon chemical 
potential (/is). The corresponding cross-over temperature range has been estimated to be 
about 170 - 190 MeV [lj]. At large QCD based model calculations indicate that the 
transition from the hadronic phase to the QGP phase is of first order. The end point of 
the first order phase transition line is the QCD Critical Point ( CP ) 

a a. 

There are 

large theoretical uncertainties of its location and even its existence is not confirmed js|. 
Experimentally, we study the QCD phase diagram by varying the colliding energy in heavy 
ion collisions j4|. The possibility of the existence of the CP has motivated our interest to 
search for it with the RHIC beam energy scan program [5!]. By tuning the collision energy 
from a center of mass energy of 200 GeV down to 5 GeV, we will be able to vary the baryon 
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chemical potential from ~ 20 to /j,b of about 500 MeV. 

The characteristic feature of a critical point are a large correlation length ( £ ) and 
critical fluctuations. Recently, theoretical calculations have shown that higher moments 
of multiplicity distributions of conserved quantities, such as net-baryon, net-charge, and 
net-strangeness, are sensitive to the correlation length £ 6]. 

In Lattice QCD calculation with /ig = 0, higher order susceptibilities of the baryon 
number, which can be related to the higher order moments of the net-baryon multiplicity 
distributions, show a non-monotonic behavior near T c [7| . A similar behavior is expected at 
CP in the finite \xb region. Experimentally, it is hard to measure the net-baryon number 
event-by-event while the net-proton number is measurable. Theoretical calculations show 
that fluctuations of the net-proton number can be used to infer the net-baryon number 
fluctuations at the CP [8j. 

In this paper we will show the centrality and energy dependence of various moments and 
moment products of net-proton multiplicity distributions from published STAR Au+Au 
data [9] at ^snn = 19.6, 62.4, 200 GeV and model calculations. 

II. OBSERVABLES 

The various moment of the event-by-event multiplicity distributions are defined as: Mean, 
M = < N >, Variance, a 2 = < (AN) 2 >, Skewness, S = < (AN) 3 >/cr 3 , and Kurtosis, 
k = < (AiV) 4 >/<t 4 — 3, where AiV = iV— < N >. Skewness and Kurtosis are used to 
characterize the asymmetry and peakness of the multiplicity distributions, respectively. For 
gaussian distributions, they are equal to zero. Thus, the Skewness and Kurtosis are ideal 
probe to demonstrate the non-Gaussian fluctuation feature as expected near the CP. In 
particular a sign change of the skewness or kurtosis may be an indication that the system 
crossed the phase boundary 

III. RESULTS 

We have calculated the various moments of net-proton (Ap = Nn~ Np) distributions 
rom transport models ( AMPT (ver.2.11) [11], Hijing (ver.1.35) [12], UrQMD (ver.2.3) 
rJ) and a thermal model (Therminator (ver 1.0) jl4j]). By using several models with 
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different physics implemented, we can study the effects of physics correlations and back- 
grounds which are present in the data and that might modify purely statistical emission 
patterns, like resonance decays, jet-production (Hijing), coalescence mechanism of particle 
production (AMPT), thermal particle production (Therminator), and hadronic rescatterring 
(AMPT,UrQMD). 

The kinetic coverage of protons and anti-protons used in our analysis is 0.4 < < 0.8 
GeV/c and \y\ < 0.5. Fig. 1 shows the number of participant (N part ) dependence of 
four moments (M, a, S, k) extracted from net-proton distributions of Au+Au collisions at 
^snn — 200 GeV for the various models. M and a show a monotonic increase with N part for 
all of the models, while S and k, which are positive, decrease monotonically with N part , which 
means the shape of the net-proton distributions become more symmetric as the centrality 
increases. The dashed lines in Fig.l are derived from the independent emission source model 
[15I ]. The centrality evolution of the various moments of net-proton distributions in Fig. 1 
can be well described by such a model. 

Fig. 2 shows the N part dependence of moment products Sg, kg 2 of net-proton distri- 
butions from the AMPT string melting model with parton cross section a pp = lOmb for 
Au+Au collisions at y/s^N = 7.7, 11.5, 19.6,39,62.4 and 200 GeV. In the upper panel, the 
Sg shows almost no centrality dependence for high energy, but it shows a small decreasing 
trend for low energies. We can also see that the Sa has a strong energy dependence decreas- 
ing with increasing energy. In the lower panel of Fig. 2, the kg 2 shows almost no centrality 
dependence and the values are around unity for all energies. 

In Fig. 3, the energy dependence of moment products 5*0", kg 2 for most central net-proton 
distributions from STAR data [9] are compared with the results from various models. We 
see the data are in good agreement with the HRG model (kb& b = 1, Sb&b = tanh(fj,B/T)) 



16| and the thermal model (Therminator) results. HIJING, UrQMD and AMPT default fail 
to describe Sa and kg 2 simultaneously. For kg 2 , the results from various models show no 
dependence on energy and are close to unity. A large deviation from constant as a function of 
N part and collision energy for kg 2 may indicate new physics, such as critical fluctuations js]. 
Recent lattice QCD calculations from [17] have shown that kg 2 non-monotonically depends 
on colliding energy in the neighbourhood of the critical point. 
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IV. SUMMARY 

Higher moments of the distributions of conserved quantities are predicted to be sensi- 
tive to the correlation length at CP and to be related to the susceptibilities computed in 
Lattice QCD and the HRG model. Various non-CP models (AMPT, Hijing, Therminator, 
UrQMD, HRG) have been applied to study the non-CP physics background effects on the 
higher moments of net-proton distributions. The moment products Sa, kg 2 of net-proton 
distributions have almost no dependence on collision centrality and Ka 2 is also found to be 
constant as a function of energy for various models. On the other hand, the high energy 
data are in good agreement with the HRG model. 

Both the model as well as the data do not show non-monotonic behavior. Thus they 
can serve as a baseline of the behavior expected from known physics effects for the RHIC 
beam energy scan. The presence of a critical point in that region may result in non-gaussian 
fluctuations and non-monotonic behavior of the observables studied here as a function of 
collision energy. 
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FIG. 1: Centrality dependence of various moments of net-proton distributions for Au+Au colli- 
sions at ^/sjvtv = 200 GeV from various models. The dashed lines represent the expectations for 
statistical emission. 
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FIG. 2: Centrality dependence of moment products Sa and no 1 of net-proton distributions for 
Au+Au collisions of various energies from AMPT String Melting model calculation. 
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FIG. 3: Energy dependence of moment products Sa and ko 2 of net-proton distributions for Au+Au 
collisions of various models and STAR data. 

FIGURE CAPTIONS 



1. Centrality dependence of various moments of net-proton distributions for Au+Au col- 
lisions at t/snn = 200 GeV from various models. The dashed lines represent the 
expectations for statistical emission. 

2. Centrality dependence of moment products Sa and na 2 of net-proton distributions for 
Au+Au collisions of various energies from AMPT String Melting model calculation. 

3. Energy dependence of moment products Sa and na 2 of net-proton distributions for 
Au+Au collisions of various models and STAR data. 



